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ABSTRACT 

The  usefulness  of  stand  density  index  (SDI)  has  been 
limited  by  lack  of  a  way  to  partition  its  nonlinear 
expression  into  additive  components  to  describe  the 
relative  stocking  of  a  stand  by  species  or  quality 
classes.     In  this  paper,  a  linear  equation  is  derived 
to  permit  such  a  partition. 

A  closely  related  expression  for  grand  fir  stocking  is 
given  that  retains  the  merits  of  SDI,  but  is  more  nearly 
proportional  to  growth  of  individual  trees  within  a 
stand.     Hence,  the  partitioning  of  stocking  into  com- 
ponents will  be  more  meaningful  to  the  forest  manager 
with  this  revised  form. 

One  important  use  of  a  measure  of  stocking  is  to  describe  the  relative 
composition  of  the  growing  stock  according  to  species,  quality,  and  size 
classes.     However,  interpretation  of  the  measured  composition  depends  on  the 
units  of  measurement.     In  a  related  paper,     I  have  proposed  that  many  aspects 
of  stocking  condition  and  control  can  be  described  in  terms  of  greater  manage- 
ment significance  when  stocking  is  measured  in  units  proportional  to  site 
utilization.     Site  utilization  can  be  measured  from  current  growth  data. 
However,  when  direct  sampling  of  growth  is  not  feasible,  the  si lviculturist 
must  fall  back  on  some  other  way  to  measure  stocking  that  approximates  the 
distribution  of  growth  among  trees  within  the  stand  of  interest.     To  be  useful 
for  describing  composition,  this  alternative  measure  of  stocking  should  be 
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additive,  tree  by  tree,  so  that  the  total  stocking  for  the  stand  will  be  the 
sum  of  its  parts.     The  purpose  of  this  Note  is  to  describe  an  additive 
stocking  measure  for  grand  fir  that  is  approximately  proportional  to  growth. 

3 

Stand  Density  Index  as  proposed  by  Reineke     is  a  useful  measure  of 
stocking  that  relates  the  logarithm  of  number  of  trees  per  acre  to  the 
logarithm  of  tree  size.     The  merits  of  this  measure  in  contrast  to  polynomial 
measures  of  stocking  such  as  basal  area,  tree-area^ratio ,  and  crown  competi- 
tion factor  have  been  demonstrated  by  David  Bruce.     He  showed  that  the 
polynomial  measures,  throughout  the  range  of  their  usefulness,  are  closely 
proportional  to  Stand  Density  Index.     Furthermore,  he  argues  that  because 
Stand  Density  Index  can  be  extrapolated  to  describe  the  effect  of  density  on 
growth  in  stands  as  young  as  nursery  seedlings,  it  must  be  the  more  generally 
valid  form  to  use  as  an  index  of  relative  stocking. 

A  disadvantage  of  Stand  Density  Index  is  that  there  has  been  no  way  to 
describe  the  contribution  of  various  classes  of  trees  in  the  stand  to  the 
total  index  for  the  stand.     The  purpose  of  this  Note  is  to  demonstrate  how 
this  calculation  can  be  accomplished,  and  to  present  a  modified  version  of 
Stand  Density  Index  for  grand  fir  (Abies  grandis   (Dougl.)  Lindl.). 


COMPUTING  STAND  DENSITY  INDEX  TREE  BY  TREE 


Stand  Density  Index  as  computed  by  Reineke  is  essentially 

/  N      2  \l/2 

log  SDI  -  log  N  +  k  log/  Z    d  /n        -  k  (1) 

\i=l    i  I 

in  wjpch  k  is  3.  constant  equal  to  1.605,  N  is  number  of  trees  per  acre,  and 
(Z  d./N)l/2  is  the  diameter  of  the  tree  of  average  basal  area.     The  index  thus 
computed  is  stated  in  terms  of  number  of  trees  per  acre  for  a  stand  of  10 
inches  average  d.b.h.     In  this  form,  the  index  is  commonly  used  to  compare 
densities  from  stand  to  stand. 

Stand  Density  Index  can  be  computed  tree  by  tree  within  a  stand.  That 
is,  the  total  SDI  can  be  subdivided  to  represent  the  contribution  to  SDI  of 
different  classes  of  trees  in  the  stand.     The  equation  for  SDI  that  is  addi- 
tive tree  by  tree  is 

N  2 

SDI  =  Z       O  +  bdm)  (2) 


or : 

N  2 

SDI  =        +  b  Z    d . 
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where  a  and  b  are  coefficients  depending  on  the  average  diameter  of  all  trees, 
and  d.  is  the  diameter  of  the         tree  in  the  tally.    The  coefficients  are 
determined  from  the  diameter  of  the  tree  of  average  basal  area  by  the  follow- 
ing equations : 

a  =  10'k  (1  -  |)    (Z  d\/N)k/2  (3) 

b  -  io-k  &  a  d./m{km-1 

Thus,  if  we  have  two  classes  of  trees  in  a  stand,  the  stocking  of  each  compo- 
nent can  be  computed  from  equation  (2)  above  by  summing  the  values  of  a  +  bd2 
for  the  trees  in  each  class.     In  turn,  the  total  Stand  Density  Index  will  be 
the  sum  of  the  partial  SDI's  in  each  class. 

A  proof  that  substituting  the  coefficients  a  and  b  as  defined  in 
equation  (3)  into  equation  (2)  yields  an  expression  equivalent  to  (1)  is 
given  in  the  appendix  to  this  Note. 

APPLICATION  TO  GRAND  FIR  SITE  UTILIZATION 

Two  modifications  of  SDI  are  here  proposed  for  application  to  grand 
fir  stands.     First,  for  stands  of  a  given  average  d.b.h.,  the  index  will  be 
rescaled  so  that  an  index  of  100  corresponds  to  the  number  of  trees  needed  to 
fully  utilize  the  site.     Second,  the  form  of  the  index  within  a  stand  will  be 
changed  from  the  form  in  equation  (2)  to  be  proportional  to  diameter  cubed. 

The  reason  for  this  latter  change  is  to  make  the  contributions  to  the 
per-acre  stocking  more  nearly  proportional  to  the  growth  of  the  individual 
trees.     In  even-aged  stands  of  grand  fir,  growth  per  tree  has  been  shown  to 
be  proportional  to  diameter  cubed.5     If  growth  is  the  criterion  for  rating 
contribution  to  stocking,  the  equation  (2)  for  SDI  given  above  overrates  the 
smaller  trees  and  underrates  the  larger  trees. 

Through  analysis  of  growth  data  from  study  plots  in  even-aged  grand  fir 
stands,5  the  following  expression  was  found  to  estimate  the  minimum  number  of 
trees  per  acre  required  to  fully  utilize  a  site  capable  of  growing  grand  fir: 

log  N  =  4.26820  -  1.76235  log^    Z^d3/N\1/5  (4) 

This  equation  is  clearly  related  to  equation  (1)  for  SDI  given  above.  The 
principal  difference  is  in  the  measure  of  average  tree  size.     The  cube  root 
of  the  average  of  the  cubed  diameters  has  been  substituted  for  the  square  root 
of  the  average  of  the  squared  diameters. 

To  emphasize  the  contribution  of    £  d  in  (4) ,  consider  this  restatement 
of  the  same  relation: 


log  I  d3  =  log  (E  d3./N)l'^k/5^  +  4.26820 

i  ^ 


(5) 


5Stage,  Albert  R.     A  study  of  growth  of  grand  fir  (Abies  grandis  (Dougl.) 
Lindl.)  in  relation  to  site  quality  and  stocking.     Ph.D.  Diss .  Univ. 
Michigan,  Ann  Arbor.  1966. 
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in  which  k  =  1.76235.     Taking  antilogs  of  (5)  leads  to 
E  d3.  -  18544.      (E  d3./N?AU5S 

Is  Is 


In  this  form,  the  right-hand  side  of  this  last  expression  gives  the 
total  sum  of  diameters  cubed  needed  (per  acre)  to  utilize  the  site  complet( 
at  the  stage  of  development  given  by  E  dj/N,  the  mean  of  cubed  diameters. 
Hence,  the  percent  utilization  of  the  site  by  trees  can  be  expressed  as 


100  E  d3 

i 


,3      o  41255 
18544.    (E  d./N)  ' 

t 


or,  simplifying: 


S  = 


E  d\ 
t 


185.44   (E  d3/N)0AUS5 

i 


Note  that  the  denominator  contains  the  average  of  cubed  diameters  over 
all  trees  in  the  stand.     Table  1  lists  the  values  of  the  denominator  for  stand 
diameters  from  3  to  30  inches. 


The  sum  in  the  numerator  can  then  be  separated  into  subtotals  over 
those  classes  of  trees  in  the  stand  for  which  stocking  fractions  are  desired. 

Obviously,  it  is  very  likely  that  stands  will  rate  over  100  percent  by 
this  method,  yet  it  is  not  possible  to  have  the  site  more  than  fully  utilized. 
The  appropriate  interpretation  of  a  stocking  index  over  100,  say,  125  percent, 
is  that  a  fraction  equal  to  25/125  or  20  percent  of  the  stocking  could  be 
removed,  yet  leave  sufficient  growing  stock  to  fully  utilize  the  site  after 
sufficient  time  has  passed  to  allow  the  remaining  trees  to  adjust  to  the  new 
level  of  stocking. 

Site  utilization  of  less  than  100  percent  could  well  be  a  desirable 
objective  in  young  stands  before  any  commercial  thinning  opportunities  exist. 
In  such  stands,  the  value  of  more  rapid  diameter  growth  might  well  exceed  the 
value  of  the  unused  site  potential. 

Example 

Site  utilization  calculations  for  a  sample  tally  from  a  fifth-acre 
plot  in  a  young  grand  fir  stand  are  illustrated  in  table  2.     In  this  example, 
the  site  is  apparently  99  percent  utilized,  about  48  percent  of  the  growth 
capacity  of  the  site  being  used  by  desirable  trees. 
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Table  1. --Values  useful  in  estimating  stocking  for  grand  fir. 


Di  amet er 

d3 

-3       0  .41255 
185.44(1  d./Ny 

3 

27 

722 

4 

64 

1,031 

5 

125 

1,359 

6 

216 

1 ,703 

7 

343 

2,061 

8 

512 

2,432 

9 

729 

2,813 

10 

1,000 

3,205 

11 

1 ,331 

3  ,606 

12 

1,728 

4,016 

13 

2,197 

4,435 

14 

2,744 

4,861 

15 

3,375 

5,294 

16 

4,096 

5,734 

17 

4,913 

6,181 

18 

5,832 

6,634 

19 

6,859 

7,093 

20 

8,000 

7,558 

21 

9,291 

8  ,029 

22 

10,648 

8,504 

23 

12, 167 

8  ,986 

24 

13,824 

9,471 

25 

15,625 

9,962 

26 

17,576 

10,458 

27 

19,683 

10,958 

28 

21,952 

11,462 

29 

24,389 

11,971 

30 

27,000 

12,484 
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APPENDIX 


PROOF  OF  TREE-BY-TREE  EXPRESSION  OF  SDI 
Equation  (2) ,  with  coefficients  a  and  b  defined  by  equation  (3)  becomes 

-V 


SDI 


=  .Z.   (a  +  bd2.)  =  Na  +  b  .Z.  d2. 
t=l  v  %J  ^=l  ^ 


io_Zc  (l  -  k/2)  [Jix  d2/Njk/2  N 


Hence 


SDI  =  10 


(i  -  w)  [   dzy/2  N  1 

1  -  K/2jj^  Ak/2 


k/ 


*  {km  cr         .z,  ^ 


10"*  (I^2  ^-fe/2 


Taking  logarithms  to  base  10 


log  SDI  =-k 


N 


+  fc/2  log  dJJ  +  (1  -  fe/2)  log  /V 


log  SDI  =  k/2  (log  .|  dt  -  log  /!/)  +  log  N  - 
log  SDI  =  log  N  +  k  log  CZ  d^-/N)1/2  -  k 


This  last  equation  is  identical  with  (1)  which  was  to  be  proven 


FOREST  SERVICE  CREED 


The  Forest  Service  of  the  U.  S.  Department  of  Agriculture  is 
dedicated  to  the  principle  of  multiple  use  management  of  the 
Nation's  forest  resources  for  sustained  yields  of  wood,  water, 
forage,  -wildlife,  and  recreation.  Through  forestry  research, 
cooperation  with  the  States  and  private  forest  owners,  and 
management  of  the  National  Forests  and  National  Grasslands, 
it  strives  —  as  directed  by  Congress  —  to  provide  increasingly 
greater  service  to  a  growing  nation. 


